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Abstract 



We consider the possibility of a successful Affleck-Dine mechanism along the 
u c d c d c direction in R-parity symmetric extensions of the minimal supersymmetric 
Standard Model (MSSM) which contain a gauge singlet superfield <fi. Such gauge 
singlets commonly occur in extensions of the MSSM, for example in models which 
seek to account for neutrino masses. We consider a two scalar Affleck-Dine mecha- 
nism, with the flat direction stabilized by a non-renormalizible superpotential term 
of the form ^^u c d c d c ~ ^0^ 3 , where if) corresponds to the gauge non-singlet flat di- 
rection. We give approximate solutions of the scalar field equations of motion which 
describe the evolution of the condensates and show that the final baryon asymme- 
try in this case is suppressed relative to that expected from the conventional single 
scalar Affleck-Dine mechanism, based on a superpotential term of the form -T^ip 4 , 



and is the supersymmetric cf) mass. It is possible for the model to generate a 
baryon asymmetry even in the limit of unbroken B-L, so long as the gauge singlet 
condensate doesn't decay until after anomalous electroweak B+L violation is out 
of equilibrium following the electroweak phase transition. This condition is gener- 
ally satisfied if all Dirac neutrino masses are less than around lOkeV. This class of 
Affleck-Dine models can, in principle, be experimentally ruled out, for example by 
the observation of a Dirac mass for the /i or r neutrino significantly larger than 
around lOkeV together with a mostly Higgsino LSP. 



by a factor 




where m s is the soft supersymmetry breaking scalar mass 



1. Introduction 

In supersymmetric (SUSY) models [lj], the occurence of flat directions in the 
renormalizible scalar potential of the minimal SUSY Standard Model (MSSM) and 
many of its extensions naturally leads to the possibility of generating the baryon 
asymmetry of the Universe via the decay of scalar field oscillations along such flat 
directions. This possibility is the well-known Affleck-Dine (A-D) mechanism for 
baryogenesis 0. Although in the limit of unbroken SUSY the renormalizible po- 
tential along these flat directions is completely flat, once soft SUSY breaking terms 
and non-renormalizible terms consistent with the symmetries of the model are added 
there will be a non-trivial potential. In the original A-D scenario it was assumed 
that the soft SUSY-breaking terms are the same as the zero temperature SUSY- 
breaking terms, which are characterized by a mass scale m s of the order of lOOGeV- 
lTeV (![]. However, it has recently become clear that the large energy density which 
exists in the early Universe will also break SUSY, resulting in soft SUSY breaking 
terms characterized by a mass scale typically of the order of the Hubble parame- 
ter H M. This large mass scale for the SUSY breaking terms radically alters the 
evolution of the scalar fields during and after inflation [|J . In the original A-D mech- 
anism, because the scalar field masses are much smaller than H during inflation, the 
classical scalar fields are overdamped and effectively frozen in at their initial values 
on horizon crossing, as generated by quantum fluctuations [0, |5|. Therefore on the 
scale of the observable Universe there is a large constant scalar field over the whole 
Universe with an essentially random phase. This then evolves into a coherently 
oscillating scalar field, corresponding to a Bose condensate with a roughly maximal 
asymmetry in the condensate particle number density. The subsequent decay of the 
condensate was shown to be easily able to account for the baryon asymmetry of the 
Universe ||. However, once mass terms of the order of H for the scalar particles 
are introduced, this picture completely changes J|. Now the classical scalar fields 
can evolve to the minimum of their potentials on a time scale of the order of H _1 . 
As a result, at the end of inflation, all the scalar fields will be at the minimum of 
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their potentials, with quantum fluctuations having an effect only on the scale of the 
horizon at the end of inflation, which is much smaller than the scale of the observ- 
able Universe. Therefore the baryon asymmetry coming from the A-D condensate 
in this case will be determined dynamically by the evolution of the scalar fields 
during the post- inflation era, with the scalar fields starting out at the minimum of 
their potentials at the end of inflation. Since the A-D mechanism is now dependent 
upon the details of the scalar potential, one has to consider each case individually 
in order to determine the magnitude of the resulting asymmetry. The asymmetry 
will be particularly sensitive to which flat direction the scalar fields oscillate along 
and to the form of the non-renormalizible superpotential terms which determine the 
minimum of the scalar potential and introduce the CP violation necessary in order 
to generate the baryon asymmetry ||]. 

In order to generate a baryon asymmetry from a condensate which decays prior 
to the electroweak phase transition (when anomalous B + L violation is in thermal 
equilibrium 0) it is necessary for the condensate to carry a non-zero B-L asymmetry. 
The lowest dimension operators which characterize the B-L violating flat directions 
in the MSSM are the dimension 2 (d=2) operator LH U and the d=3 operators u c d c d c , 
d c QL and e c LL [|J. (These operators may be thought of as the superpotential terms 
responsible for lifting the flat directions or as the scalar field operators which are 
responsible for introducing explicit B-L violation into the scalar field equations of 
motion. These are naturally connected by the relationship between the soft SUSY 
breaking terms and the superpotential terms JIL M). These operators characterize 
the flat directions in the sense that the scalar field operator characterizing a par- 
ticular flat direction will have a non-zero expectation value along that direction. 
We will refer to the flat direction which gives a non-zero expectation value to LH U 
and to u c d c d c as the "LH U direction" and the "u°d c d c direction" respectively. (The 
A-D mechanism along the d c QL and e°LL directions will be essentially the same as 
that along the u c d°d c direction; we will concentrate on the u c d°d c direction in the 
following). The LH U direction and the u c d c d c direction are orthogonal in the sense 
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that they cannot both be flat simultaneously The LH U direction has recently 
been considered by a number of authors 0, 0]. In the present paper we will focus 
on the u c d c d c direction. 

The simplest implementation of the A-D mechanism along the u c d c d° direction 
would involve adding to the MSSM superpotential a d=3 term of the form u°d c d c . 
However, this term would be phenomenologically dangerous, as it would introduce 
large B violation into the MSSM unless its coupling was extremely small. (For 
a review of the constraints on B and L violating couplings in the MSSM see ref- 
erence For example, squark mediated proton decay imposes the constraint 
| A' A" | < 10 -24 for the light quark generations, where A' is the d c QL coupling and 
A" is the u c d°d c coupling ||. Such dangerous B and L violating terms are usually 
eliminated from the MSSM by imposing R-parity (R p ) [p], [8]. Imposing R p implies 
that the first B-L violating operator in the MSSM which is nonzero along the u c d°d c 
direction is a dimension 6 operator, u c u°d c d c d c d c ||. However, as discussed in ref- 
erence [|J (and briefly reviewed in the present paper), in A-D models where the 
natural scale of the non-renormalizible terms is the Planck scale, d = 6 A-D models 
can have an acceptably low B asymmetry only for very low reheating temperatures 
Tr, #Tr < lOGeV, where 9 is a CP violating phase, d = 4 A-D models, on the 
other hand, can be compatible with a much wider range of reheating temperatures, 
up to 10 9 GeV or more 0]. Thus with only the particle content of the MSSM, the 
u c d°d c direction would be disfavoured in the simplest models (those based on Planck 
scale non-renormalizible terms) relative to the LH U direction, which can utilize the 
Rp-conserving d=4 operator (LH U ) 2 . However, if we were to consider extensions 
of the MSSM which involve the addition of an R p -odd gauge singlet superfield <fi, 
then we could form the R p conserving d=4 operator 0u c d c d c . The addition of such 
a gauge singlet superfield to the MSSM is a very common and natural feature of 
many extensions of the MSSM. In particular, in models which seek to account for 
neutrino masses, the gauge singlet superfield would correspond to a right-handed 
neutrino superfield. It is the purpose of the present paper to determine whether it 
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is possible to generate the observed B asymmetry along the u c d c d c direction via the 
operator 0u c d c d c and, if so, to compare the resulting asymmetry with that coming 
from the more conventional LH U direction. 

The paper is organized as follows. In section 2 we discuss the model and the min- 
imization of its scalar potential. In section 3 we consider the scalar field equations 
of motion and the formation of the coherently oscillating scalar field condensates. In 
section 4 we discuss the condensate particle asymmetries and the resulting baryon 
asymmetry. In section 5 we discuss the constraints on the reheating temperature 
after inflation. In section 6 we discuss the thermalization and decay of the conden- 
sates and the upper limits on Dirac neutrino masses in the limit of unbroken B-L. 
In section 7 we give our conclusions. 
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2. d=4 Affleck-Dine mechanism along the u c d c d c direction 
We will consider throughout the simplest scenario, in which it is assumed that 
inflation occurs with an energy density consistent with the density perturbations 
observed by COBE, corresponding to H ps 10 14 GeV ||, with the inflaton $ subse- 
quently undergoing coherent oscillations about the minimum of its potential. We 
will also require that the reheating temperature Tr is low enough not to thermally 
regenerate gravitinos [[HJ], which implies that Tr is less than about 10 10 GeV, corre- 



sponding to H not much larger than lGeV. After reheating we will assume that the 
Universe is radiation dominated throughout, with no further significant increase in 
entropy. In general, when inflation ends H will not be much smaller that its value 
during inflation (even in $ 2 chaotic inflation, the value of H when the inflaton $ 
starts oscillating is not much smaller that 10 13 GeV |11|). Therefore the coherent 
oscillations of the A-D field, which will begin once H ps m s ps lOOGeV, will begin 
during a matter dominated era, with the energy density of the Universe dominated 
by inflaton oscillations [|J. 

It is now understood that in most supergravity models, the energy density that 
exists in the early Universe will break SUSY, introducing soft SUSY breaking terms 
characterized by a mass scale typically of the order of H ||. We will therefore 
consider in the following soft SUSY breaking terms of the form 

V soft = (m 2 - CiH 2 )!^! 2 + (B a W 2 a + h.c.) + (A a W n a + h.c.) (2.1), 

where W 2 a are superpotential terms bilinear in the fields and W n a are terms of order 
n in the fields. A a and B a are defined by A a = A as + a a H and B a = B as + b a H, where 
A a s and B a s may be thought of as the zero temperature soft SUSY breaking terms 
from a hidden sector of N=l supergravity whilst a a H and b a H are due to SUSY 
breaking by the energy density in the early Universe. We will assume throughout 
that A as ps B as ps m s . We will also assume that a 2 ~ b 2 ps |c a |. In most supergravity 
models we expect that |c a | ps 1 [[|, although in some models |c a | may be smaller; 
for example |c a | ps 10~ 2 occurs in supergravity models with a Heisenberg symmetry 
[p~2]j . Then for the case with c a > the A-D scalar, corresponding to a renormalizible 



flat direction in the scalar potential, will have a non-zero value at the minimum of 
its scalar potential at the end of inflation, with the potential being stabilized by the 
contribution from the non-renormalizible terms in the superpotential H|. 

To implement the d=4 A-D mechanism along the u c d c d° dirrection we will 
consider an R p symmetric extension of the MSSM defined by the superpotential 
W = W sm + W' + W„, where W sm is the MSSM superpotential, W' is defined by 

nu , A 



W 



+ -0uW 



V 
4M' 



(2.2) 



and W„ is given by 



A„0H u L (2.3). 



In addition to these terms we would expect terms of the form 0d c QL and 0e c LL. 
For simplicity we will not include these terms explicitly. The operator u c d c d c = 
e a/ 3 7 u^d^d^ (where a, [3 and 7 are colour indices and generation indices are, for 
now, suppressed) is antisymmetric in the d c scalar fields. Therefore the d c should 
be from different generations, which we will denote as d c and d° '. With u c , d c and 
d c ' having different colour indices, the F-term contribution to the scalar potential is 
then given by 



Vf = E 



aw 



m!|0| 2 + 



IAI 



IAI 



M 2 



M 2 

|d c d c 'l 2 + |u c d c | 2 



u c d c d c '| 2 



\V\ 
M 2 



I0| f 



|u c d c '| 2 



+ 



m 



'^llW' + ^) + ^(uW')^ + h c 



(2.4). 



The direction with only < u c >, < d c >, < d c ' > and < <fi > non-zero is F-flat in 
the MSSM, with only the terms in W' lifting this flatness. The D-term contribution 
to the scalar potential, 

V D = E|"l $tT r $ | 2 (2-5), 
where the $ are the multiplets of the gauge group i with generators Tf , also vanishes 



so long as u c , d c and d c have different colour indices and |u 



c|2 



ic|2 



ic'|2 



The phases of u°, d° and d° are not, however, fixed by the MSSM F- and D-flatness 



conditions. We will see in the following that the important phases for the A-D 
mechanism are 5 V and 84, where < u c d c d c ' >= v 3 e 1<5v and < >= v^e'V 
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2.1. Potential minimization in the 77 — > limit 

We first note that in the SUSY limit with H=0 and with 7^ 0, there is a 
minimum with v^ 7^ which is degenerate with the = minimum and which 
could be phenomenologically dangerous. To be precise, from equation (2.4) we see 
that the SUSY minima correspond to 4> — and (p 2 = — (with u c = d c = d c ' = 
at both minima). The <fr 7^ minimum results in a dangerous u c d c d c superpotential 

|AI / |m I \ 

term with coupling y^^J • (In general we would expect similar couplings for 
the d c QL and e c LL terms). Thus, if we consider M < Mpi (where Mpi is the Planck 
scale) and |m^| > m s > 10 2 GeV, then we see that (^-f^ > ^lO^ 8 , which, 
for |A| and \rj\ not much smaller than 1, would result in an unacceptable squark 
mediated proton decay rate ||. In order to avoid this danger we must therefore 
ensure that there exists a minimum of the scalar potential which has v^O for large 
H but which evolves to the v = v^ = minimum as H tends to zero. 

With regard to the scale of the non-renormalizible terms M, we will set this 
to equal the Planck scale by convention. Then the coupling A can take values 
small or large compared with 1, depending on the natural mass scale of the non- 
renormalizible terms. Values much larger than 1 would correspond to the case 
where the natural mass scale of the non-renormalizible terms is much smaller than 
the Planck scale, for example a grand unification scale. On the other hand, if the 
natural mass scale of the non-renormalizible terms was of the order of the Planck 
scale, then we would expect that |A| < 1. 

Typically, we would not expect A and r\ to be very different in magnitude. How- 
ever, we would like to be able to minimize the potential analytically. We find that 
we can do this for the case of 1 77 1 small compared with |A| {\r\\ < 0.1|A| is sufficient), 
in which case it may be shown that the terms in the scalar potential proportional 
to 77 can, to a good approximation, be neglected. We expect that the possibly more 
likely case with \v\ « |A| will be qualitatively similar. In the following we will 
consider the minimization of the potential in the 77 — * limit. 
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In the ?7 — > limit the scalar potential becomes 



,|AP 



IAI 



V = (K| 2 + m 2 - c,H 2 )v 2 + K - c v H 2 )v 2 + 3^v 2 v 4 + ^v b + 



M 2 



+ 



m^A 
M 



v^vV^ - ^ + h.c. 



+ 



M 



v v 3 e i(,5v+ ^ ) + h.c. 



5^v!e 2 ^ + h.c. 



(2.6). 



By an choice of the phases of the scalar fields we can make and A real. We may 
also choose m^ to be positive. To a reasonable approximation we can neglect the 
term proportional to B^,. This only contributes a term of the order of m0(m s +a i H)v 2 6 , 
which is less than or of order of the (m 2 . + m 2 — C0H 2 )v 2 term. The phases 5$ and 5 V 
will then adjust to minimize the cross-terms (where we use "cross-terms" to denote 
terms which are the sum of a term and its hermitian conjugate). The potential will 
then have the form 

V = (m 2 + m 2 - c^H 2 )v 2 + (m 2 - c v H 2 )v 2 + — v 2 v 4 + — v 6 - 2m^ V<p v 3 (2.7), 
where fh^ = m^ + jA^j. In general, the minimum of this potential is given by 



v 



IAI 



m v° 



M(( m | + m 2- c ,H 2 ) + fgv 4 ) 



(21 



and 



v(v 4 - ^v,v + 2vjv 2 + ^(m 2 - c v H 2 )) = (2.9). 



We next consider how the minimum evolves from an initially large value of c^H 2 . 
(i) c H 2 > m 2 + m 2 



In general the minimum in this case is given by 

a* M 2 H 2 



(2.10) 



and 



(1-a) 3A 2 
1 rrv(l-a) 1/2 v 1 a A (1 - a) 1 ' 2 



9 V3E cJ/2 
where a is the solution of 



a 



a 



(2.11), 



a 3 -(l----)a 2 -|-« + i-= (2.12). 



2 a 2 
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as may be seen by taking the H 2 terms large compared with the mass terms and 
substituting (2.10) into (2.9) and (2.8). Typically (1 — a) is of the order of 1. For 
example, if a 2 is small compared with c v then a = (1 — ^-), whilst if a 2 = c^ = c v 
then a = —2. Therefore we can roughly say, with c v « w a 2 , that 



v 



i/4 (MH) 



1/2 



IAIV2 



(2.13). 



Thus v and are initially of the same order of magnitude, 
(ii) c^H 2 < m 2 + m 2 

In this case we find that it is consistent to assume that m 2 . 
which case is given by 



m 



2 > 3Af_. 

M 2 



m 



IAI 



m^v" 



M ( m 2 + m 2) 

Solving (2.9) for v, we find solutions v± given by 



(2.14). 



2\2 



V± = "4 A 2 



1 M 2 (m 2 + m 2 ) 



1 



m 2 



(m 2 +m s 2 ) 



1 ± 



1 - - 



~ 2 I 2 

ml m 



c v H 2 ) 



3 (m 2 + m s 2 - m 2 ) 2 



1/2' 



(2.15). 



For c v H 2 > m 2 , vi is negative, and so there is a minimum at v + with no barrier 
between v = and v = v + . Once c v H 2 is less than m 2 , a barrier appears at v_. The 
minimum at v + subsequently becomes unstable once 



3(2|A A K + |A A | 2 ~™ 

m — c v ri > 



m: 
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(m^ + |A A |) : 



(2.16). 



For m^ large compared with |A As | and m s we see that we must have m 2 > ||Aa s | 2 
in order that the dangerous v + ^ minimum becomes destabilized as H — > 0. As 
m^ — > 0, this condition becomes m 2 > | TaTT^ • Thus we see that it is non- 
trivial for the dangerous v^O minimum to become destabilized as H — > 0. So long 
as the potential can be destabilized, however, it will generally destabilize once m 2 
is greater than c v H 2 up to a factor of order 1. The values of the fields when the v + 
minimum becomes unstable are then given by 



1M 2 
4 A^ 



'm|+m^ 
ml 



^(mj-mj-m 2 ) (2.17) 



11 



and 

l (m 2 -m 2 -m s 2 )V 2 

Noting that |A>| s;m s for c v H 2 < m 2 , we find that, for > m s , 

1 /m \ 1/2 

v ^7lU) v (2 ' 19) - 

whilst for < m s m v. So for > m s we find that becomes suppressed 
relative to v. 
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3. Bose condensate formation 

To discuss the formation of the Bose condensates, that is to say, the way in which 
the scalar fields start oscillating freely about the minimum of their potentials, we 
consider the equations of motion of the scalar fields. Strictly speaking we have four 
scalar fields; u c , d c , d c ' and <p. However, since we are considering the evolution of the 
classical fields along a D-flat direction, we may impose that the u c , d c and d c ' fields 
have the same magnitude. Although the phases of these fields could be different, 
the equations for u c , d c and d c ' are identical and the initial values of the fields are 
the same. Therefore we may assume that their phases remain equal throughout. 
The equations of motion are then given by, 

i, K + 3H^r = -[(m s 2 - c v H 2 )^ R + w <fa& + ^ R 

-/9^r0r - a(t)(c e R V4 + s e (> R 0i + 2^r0r^i))] (3.1), 
A + 3H^ = -[(m 2 - c v H 2 )V>i + W <Pl^i + W2 KA 
-Ptylfa - 2V>rV>i0r) - a(t)(so<M& ~ cei^i + 2^ R R ^i))] (3.2), 
R + 3H0 R = -[(m 2 + m s 2 - c^H 2 )0 R - « + JpK^R 
-7(c e R + s e 0i) - a(c e V R + Sse^R^i))] (3-3), 

and 

3A 2 

fa + 3H0j = -[(m 2 + m s 2 - c^H 2 )^ - 3/3^1 + y^V^i 
-7(s e R - c e 0i) - a(s^R - Scgip^ipi))] (3-4), 
where tp represents the u c , d c and d c ' fields and we define a, (3 and 7 by a 



M 



(3 = and 7 = iB^m^l, where — ^p- = —(3 (we choose this to be real and negative 
by an choice of phase), (7^7)^ = —ae 10 and (B^m^ = — 7e le . c e (c e ) and sg (s e ) 
denote Cos 9 (Cos e) and Sin 9 (Sin e) respectively. In writing these equations we 
have assumed that the real parts of the fields are large compared with the imaginary 
parts, which turns out to be a reasonable approximation. Throughout our discussion 
of the equations of motion we will focus on the most likely form for the soft SUSY 
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breaking terms in the early Universe, corresponding to the 
II- 

Before the v ^ minimum is destabilized, the fields will be at the minimum of 
their potentials, corresponding to setting the right hand side (RHS) of the equations 
of motion to zero. The v ^ minimum at v = v + will become destabilized once 
c v H 2 < m;? and the tfin field will begin to roll once H 2 < m 2 . Once ^r starts to roll, 
we will see that the other fields follow the minimum of their potentials as a function 
of ^R,(t) and oscillate about this minimum until they become freely oscillating about 
the = minimum of their potentials. 

We first consider the evolution of the real parts of the fields, beginning with 0r. 
We consider the solution of the 0r equation of motion in the limit where the terms 
proportional to V'r^R) 7 an d are neglected, which turns out to be a reasonable 
approximation. We will also treat 9 as a time-independent constant. (We will 
comment on this later). The </>r equation of motion is then approximately given by 

R + 3H0r « -[(m 2 + m 2 - c H 2 )0 R - {ac e + (3.5). 

From now on we will neglect the C;H 2 mass terms in the equation of motion as these 
will quickly become negligible as the Universe expands. We will also neglect the 3H0 
and 3H-0 damping terms, since we are considering m s > H. The effects of damping 
will, however, be included in the time dependence of the amplitude of oscillation of 
the fields. 

We wish to show that in the solution of this equation </>r oscillates around the 
minimum of its potential as a function of ifa. To see this let 0r = R + 50r, where 

0R=^4^-#i (3-6) 

m <f> T U1 S 

is the minimum of the potential as a function of ifet(t). Then the £0r equation of 
motion is given by 

+ 6#rV> 2 ~ -(m 2 + m 2 )50 R + 3^m s 2 ^ (3.7). 

In this we have used -0r ~ — m 2, 0R, which will be shown later to be true. 50r 
will then grow from an initial value of zero until the mass term on the RHS of 
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(3.7) proportional to <50 R becomes dominant, after which 5</> R will oscillate about 
the minimum R = </> R with frequency ~ (m| + m 2 ) 1 / 2 , the terms on the RHS 
proportional to tp^ being rapidly damped by the expansion of the Universe. The 
initial value of the £0 R oscillation amplitude is therefore given by <50 Ro , where 

( 3 - 8 )- 

+ m| 

It is straightforward to show that ~ m = which is less than or about equal to 

to <f>no m^+ m 2' H 

1. Eventually the amplitude of &/> R will become larger than that of R , in which case 
5(/) R p» R will effectively oscillate freely around R = with frequency (rn^ + m 2 ) 1 / 2 . 

We next consider the solution of the ipn equation of motion. On introducing 
4>r = + R , we find that, for m 2 > c v H 2 , if) R begins oscillating with a frequency 
approximately equal to m s . This is not at first obvious since the V'r^r an d 
terms on the RHS of the -0 R equation of motion are initially large (~ (m s m^)^) 
compared with the m 2 -0 R term. However, it turns out that there is a cancellation 
between these higher-order terms, such that the sum of these terms on the RHS of 
(3.1) contributes initially only ~ m 2 ^ R and then rapidly becomes small compared 
with the m 2 -0 R term as ^r decreases with the expansion of the Universe. Therefore 
ipn will essentially oscillate with frequency approximately equal to m s once the 
v = v + minimum becomes unstable. 

Thus we can summarize the evolution of the real parts of the fields by 

Vfc(t) w A v ,(t)Cos(m s t) (3.9) 

and 

</> R (t) = R + 5<j) R « #R(t) + A (t)Sin ((m + m s )t) (3.10), 

where the time dependence of A^(t) and A^(t) is due to the expansion of the Uni- 
verse during inflaton matter domination, A(t) oc a(t)~ 3//2 , where a(t) is the scale 
factor. 

We next consider the evolution of the imaginary parts of the fields. The evolution 
of 0i and ip\ is similar to the evolution of R i.e. they follow the minimum of their 
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potentials as a function of ipR(t). We first consider the solution of the ipi equation 
of motion. 

In the ipi equation of motion we may roughly absorb the terms proportional to 
acg into the terms proportional to (3. Thus we may neglect these terms for now. 
We will also set the Pip^i term to zero for now. With these assumptions only the 
phase 8 contributes to the imaginary parts of the fields. We will comment on these 
assumptions later. 

Suppose the i/j-r field starts oscillating at t G . (For convenience we will set t G equal 
to throughout). Initially by a choice of the phase of the scalar fields, we can set 
the phase 8 to zero at t Q . The subsequent evolution of the phase 8(t) is then found 
from 

ae .,.) = _( (A.. + a, He^)A y (1+a>oe , rl ^ 

where a is the phase difference between the A A s and a a = a a e 1CT terms. Since 
during matter domination 

H = - ^ (3.12) 

(l + fH t) 1 ; 

where H Q = H(t ) ~ m s , we see that the phase 8{t) will reach its maximum roughly 
during the first oscillation cycle, in a time St 1 ps m^ 1 . Since the conden- 
sates will form during the first few oscillations of ifa, it is a reasonable approximation 
to set 8(t) to its constant maximum value throughout. With the above assumptions 
the ipi equation of motion can be reasonably approximated by 



A 2 



(3.13). 



From this we see that the minimum of the ipi potential as a function of ^r is given 
by 

-j- as^M 2 

A2 ( l + W) 

Thus, noting that, at H ~ m s , ™|^ 4 is small compared with 1, we see that ipi 
is initially proportional to ifa. This will continue until ^r decreases during its 
oscillation to the point where > 1. 
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To understand the evolution of fa let fa — fa + Sfa. Substituting into the fa 



equation of motion, we find that, for < 1, Sfa satisfies 

Sfa « A2 m s ^R - (3-15). 

Thus will grow from Sfa = to a value given by 

as e ?7m 2 M 4 

(3.16). 



We see that the condition < 1 is equivalent to ^ ^i- Initially ch/>i has a 

value 

sg|A A |(|A A | +m^) / m s \ 

o^io- 7—2- ^ : Ufc (3.17). 

(m^ + m 2 ) ^m^ + ms/ 

Thus initially ~ ( m "| s m ) s#, which is small compared with 1 for m^ large 
compared with m s or 9 small compared with 1. During the subsequent -0r oscillation, 
we see that, so long as < 1, <fr will be proportional to ^r. Therefore fa will 

initially be in phase with However, as V'r decreases, for a period St during 
the fa oscillation 5 fa will become larger than fa and the approximate equation of 
motion (3.15) will no longer be valid. During this time the m 2 ^ term in the fa 
equation of motion will dominate and the fa oscillation will continue with frequency 
m s . However, since during this period the effective mass term in the fa equation 
of motion will differ from m s by a factor of order 1, it will be possible for the fa 
phase to shift relative to fa& by approximately m s 5t. As a result, for a fraction m s St 
of the total fa oscillation, there will be a phase shift approximately given by m s St. 
Thus the average phase shift between fa and V>r over the period of the fa oscillation, 
which, as discussed in the next section, is relevant for determining the ip asymmetry, 
is given by S p ps (m s 5t) 2 . St corresponds to the time during which -0 R < . For 
a matter dominated Universe, with ip^ oc a(t)~ 3 / 2 , we find that S p ~ (m s 5t) 2 ~ 

/ \l/2 / \2 / \ 1/4 

V m"+mJ (it) ■ S p caches its largest value, S p ps 1, once H ps ( m "j s m J m s . 
We also note that for 9 small compared with 1, ps ^p- ps ( — ^ — ) Sg is small 
compared to 1, as has been assumed throughout. 

These results hold if (i) acg is small compared with (3 and (ii) if the Pfaijj^ term 
in the fa equation of motion can be neglected. The effect of the terms proportional 
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to acg will be to effectively multiply the terms proportional to (3 by a factor ~ 
(1 + — eg). However, we can still reasonably ignore the /3-0R-0I0R term in the fa 
equation of motion, even with this factor. The effect of including the fifaip^ term 
in the fa equation of motion (together with the factor from (i)) turns out to be to 



approximately replace Sg by s^, where = sg + (^ m J s e . Thus an imaginary part 
for ip can also be generated by the phase e so long as has a mass term in the 
superpotential. 

We finally consider fa. With fa = fa + 5 fa, where 



ml + m 



S / 



and 

3 (A A c e + 
(m s 2 + m 2 

the fa equation of motion can be written approximately as 



« kasg + 777s. - 2 2 2 , 2 kas e + 777s. 3 
501 + 6 j^V'rV'r ~ - m + m s 50i + 3m — R 3.20 . 

( m + m s) ( m + m s) 

Thus 50i will increase from zero to an initial oscillation amplitude given by 

(m 2 +m s 2 ) 2 

and will subsequently oscillate freely with frequency (m 2 + m 2 ) 1 / 2 . In general 
50i o ^ 0i o- It is straightforward to show that 50r and 50i reach their maximum 
values and begin oscillating in a time 5t ~ (m^ + H1 2 )" 1 / 2 < m" 1 and so will become 
freely oscillating within the first few oscillations of the -0r field. 

It is important to emphasize that it is the oscillations of 50r and 50i about the 
R and fa minima which evolve into the Bose condensate as R and fa become 
smaller than 50r and 50i. This occurs once H < ( m ™ s m ) m s- We will show in the 
next section that the average asymmetry over the ip oscillation period m s comes 
purely from 50 and not from 0. 

In summary, the if) and fields will begin to oscillate at H ps m s , with the 
initial values of the oscillating scalar field amplitudes relevant for the formation of 
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condensate particle asymmetries given by 

M 1/2 

feo ~v+ « p^(m s 2 + m s m^/ 4 (3.22), 

/ 377m 2 \ „ m 2 (m 2 + num s ) 1/2 

2 2 M^Ro- V ; fee (3.24), 



ml + m 2 y (m^ + nx 



s; 3 



and 



(kas e + 7 7ys e ) 3 ^ Sg-m^K 2 + m^m s ) 1 / 2 
^ Io ~ K + m 2 ) ^° (m 2 + m 2 ) 2 (3 ' 25) ' 



where in the final expressions we have set |A A | |B^,| ~ m s in order to show the 
dependence on m s and m^. However, as noted above and discussed further in the 
next section, although the ip scalar begins oscillating at HKm S) for the case where 
m^ > m s the full phase difference 5 P between and ifa and the associated ift particle 

/ \l/4 

asymmetry will only form once H is smaller than ^ m< ^f m J ni s . 

The approximations used in obtaining these results are good for ^> m s . For 
m^ < m s some of the assumptions made are only marginally satisfied or even slightly 
violated (although not strongly violated). However, we expect that the above results 
for the initial amplitudes will still be qualitatively correct, giving the correct order 
of magnitude for the resulting particle asymmetries. 
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4. Particle and Baryon Asymmetries. 

In the limit where we retain only the mass terms in the ifj and <fi equations of 
motion, there is a global U(l) symmetry. This is broken by the B and L violating 
terms coming from the non-renormalizible terms in the scalar potential, which give 
rise to a non-zero U(l) charge in the condensate, corresponding to an asymmetry 
in the number of ip and (f> particles. We first consider the ip asymmetry. The 
asymmetry in the number density of ip particles is given by 

= i(V* - ^ VO = -2(VW>i - V#r) (4-1). 

Note that in the case where ?/>r and ipi oscillate with the same frequency, there 
must be a phase difference between the oscillating t/'r and ipi fields in order to have 
a non-zero asymmetry. As discussed in the last section, there is a time-dependent 
phase difference 5 P between -0r and ipi (averaged over the period of oscillation m^ 1 ), 
which for the case S p < 1 is given by 




and by 5 P ~ 1 otherwise. We can characterize the ip particle asymmetry by the 
asymmetry at H « m s that would evolve to the correct ip asymmetry at present, 

ity w 28 p m^ K oi>\ o (4.3) , 

where the value of 5 P is determined by the value of H at which the condensate 
thermalizes or decays. 

For the case of the (f> asymmetry, since, for m^ > m s , the 50r,i oscillate with a 
greater frequency than the (f>^i, substituting <p = <fi + 5<p into (4.1) shows that the 
non-zero average asymmetry over the period H _1 > m^ 1 > (m 2 + m 2 ) -1 / 2 will be 
that purely due to <50r and 54>\- Thus we find that the number density asymmetry 
in particles is given by 

n</> (m + m s ) S p 
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The suppression of relative to n^, when > m s is significant, since in the limit 
where — > and S p — > 1 we may define an unbroken B-L asymmetry by assigning 
B=l to <j). This is broken by m^, suppressing n^ and so preventing any possibility 
of a cancellation between the B-L asymmetries coming from decay of the ip and 
condensates. 

We next compare the ip asymmetry from the above ^0^ 3 -type model with that 
expected from the more conventional A-D mechanism based on a single field with 
a B violating term of the form -^^> 4 together with the above form of H corrections 
to the SUSY breaking terms. In this case the initial values of the ipn and ipi fields 
are given by 



/ 1 x 1 / 4 KM) 1 / 2 , x 

(4 - 5) ' 



and 



^io^ (4.6), 



m 



2 



with tp R and tpi subsequently oscillating with a phase difference of the order of 1. 
Thus in this case we find that 

with the asymmetry being fully formed at H ps m s . Therefore, for the case where the 
CP violating phases sg and s^ and the coupling A have the same values in both cases, 



we find that the asymmetry from the ^0-0 3 -type superpotential term is related to 



A x,/,3 

aoy iiiiiiciji v iiuiii inc 

that from the term by 

/ \ 1/2 

n V-~^ P ~ X~ n ^o (4.8), 
\ m^ + m s J 

where n^, is the asymmetry expected from the conventional single field A-D mecha- 
nism. Thus we see that, even with 5 P ~ 1, for m^ > m s there is a suppression of the 

/ \ 1/2 

ip asymmetry by a factor (^ m ^°j! m J relative to that expected from the conventional 
single-field A-D mechanism. 

The baryon asymmetry from the tp condensate is simply the asymmetry n^, 
multiplied by a factor for the number of baryons produced per ip decay. Since we 
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have in fact three condensates, corresponding to u c , d c and d c ', each of which carries 
baryon number 1/3, we see that ne = ity. The baryon asymmetry to entropy ratio 
after reheating is then found by noting that the ratio of the if) number asymmetry 
to the inflaton energy density during inflaton oscillation domination is constant. 
Since the reheating temperature is given by p\ = kpT^ and the entropy density by 
s = k s T| (with k T = ^§fl and k s = ^), where g(T) = g b (T) + |g f (T) and 
gb(T)(gf(T)) are the number of bosonic (fermionic) degrees of freedom in thermal 
equilibrium at temperature T ||), it follows that the baryon-to-entropy ratio is 
given by 

=2 = =-^ (4.9) 
s k T T R pi 

Thus, with the energy density dominated by inflaton oscillations when the Bose 
condensate forms at H ps m s and with ne = m s e s p™ sM ^_JHs__^ ^ ( we find that 



n B Sit SqO p / m s \ 



pi 3 |A|Mpi \m s + m^ J 
where we have used M = Mpi and |Aa| ~ m s . Therefore 

n B 2irsn5 p ( m s \ 1/2 T R 



(4.10), 



(4.11). 



s |A| \m0 + m s y Mpi 
Comparing this with the observed asymmetry, — m 10 -10 , and noting that the 



thermal gravitino regeneration constraint implies that ^ 10~ PH , we see that 
m<£ cannot be too large compared with m s if is not large compared with 1, as 
we would expect if the natural scale of the non-renormalizible terms was less than 

s-S 

or of the order of M P i. On the other hand, even if -tp is large compared with 1, 
the suppression factor [ m ^ m J would still allow Tr to be consistent with the 
observed B asymmetry for values of Tr up to the gravitino constraint, so long as 
was sufficiently large. 

Thus we can conclude that, in the R p symmetric MSSM extended by the addition 
of a gauge singlet scalar, it is indeed possible to have a successful d=4 two-scalar 
0^ 3 -type A-D mechanism along the u c d c d c direction. The resulting baryon asymme- 
try, assuming that the asymmetries are able to fully form and do not thermalize or 
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decay before the phase S p (t) has reached its maximum value (we discuss this possi- 
bility in the next section), receives an overall suppression by a factor approximately 

/ \ 1/2 

V m"+m J relative to that expected in the case of a conventional d=4 -0 -type A-D 
mechanism based on a single A-D scalar field. This suppression factor allows for a 
wider range of non-renormalizible couplings and reheating temperatures to be com- 
patible with the observed baryon asymmetry than in the case of the conventional 
A-D mechanism. 
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5. No-Evaporation Constraint 

An important constraint on the reheating temperature comes from the require- 
ment that the interaction of the condensate scalars with the radiation energy density 
due to inflaton decays prior to reheating does not lead to the condensate thermaliz- 
ing via scattering with the background plasma before the particle asymmetries can 
be established ||J]. We refer to this as the no-evaporation constraint. 

We first consider the ip condensate. In general, there are two ways in which 
the condensate can be destroyed at a given value of H; thermalization and decay. 
Thermalization of the condensate will occur if (i) the rate of scattering of the thermal 
plasma particles from the condensate scalars, T s , is greater than H and (ii) if the 
mass of the initial and final state particles are such that the scattering process is 
kinematically allowed and the scattering particles in the plasma are not Boltzmann 
suppressed. Decay of the condensate via tree level two-body decays will occur if 
(i) the decay rate of the condensate scalars in their rest frame, Td, is greater than 
H and (ii) if the final state particles, of mass < ip >, where A^, is a gauge or 
Yukawa coupling and < if) > is the amplitude of the ip oscillation, are lighter than 
the condensate scalars. Since the time over which the real and imaginary parts of 
the fields start rolling and so the asymmetries start to develop is of the order of H _1 
at H m s , we must ensure that thermalization and decay does not occur on a time 
scale small compared with m" 1 . 

Prior to reheating, the radiation energy density coming from inflaton decays 
during inflaton matter domination corresponds to a background plasma of particles 
at a temperature T r , where, assuming that the decay products thermalize, T r is 
given by @, § 

T r « kr(MpiHT^) 1 / 4 (5.1), 

where k r = (^l^) 1 ^ « 0.4 (using g(T) « 100). 

For the case of the i/j condensate, the particle asymmetry begins to form at 

/ \ 1/2 / \ 2 / \ 1/2 

H « m s , at which time o p (t) ~ ( m + S m</1 J \Ti) ~ ( m +m^, J ' an< ^ subsequently 
grows until o D ~ 1 at H ~ ( — — m s . 
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Thermalization of the ip condensate is possible if the rate of scattering of the 
plasma particles in equilibrium at temperature T r from the condensate scalars is 
sufficiently large. For the case of t-channel scattering of condensate scalars from 
plasma fermions via SU(3) gauge boson or Yukawa fermion exchange interactions, 
the scattering rate is given by T s as krcrA^T r , where is the gauge or Yukawa 
coupling, a = x ^ 2x ) f° r ^ ne g au S e boson exchange and a = |log (j^j for the Yukawa 

2 

fermion exchange, with x = ^ and hia the mass of the exchanged gauge boson or 
fermion, and where for scattering from a single Dirac fermion kr ~ ps 3xlCT 3 . 
The condition for the plasma to be able to thermalize the condensate is then that 
T s > H, which implies that 



We see that this will be satisfied by the gauge couplings, the top quark Yukawa 
coupling and, marginally, the bottom quark Yukawa coupling. Thus, in order to 
prevent the thermalization of the condensate, we must require that < ip > 
> 3T for these couplings, to ensure that the associated scattering processes are 
kinematically suppressed. To be safe, we will conservatively require that 
X^p < if) > > 30T, in order to suppress scattering from plasma particles with energy 
larger than the mean thermal energy 3T. For the case of the u c d c d c direction, we 
note that the smallest Yukawa coupling to the condensate scalar if), which will 
typically involve a linear combination of all three generations of down squark, will 
equal the b quark Yukawa coupling up to a factor of the order of 1. Thus, in 
general, kinematically suppressing the b quark Yukawa interaction will ensure that 
the if) condensate is not thermalized. Although the b quark Yukawa coupling, Ab ~ 
4xl0~ 2 for the MSSM with equal Higgs doublet expectation values, only marginally 
satisfies the condition for thermalization, and so might allow the condensate to 
form without large suppression of the asymmetry, if we were to consider t-channel 
fermion exchange scattering from light gauginos and squarks in the plasma, then the 
constraint (5.2) would apply to the combination A^ ~ (gAb) 1//2 , where g is a gauge 
coupling. This would exceed the lower bound (5.2). Thus we will conservatively 




(5.2). 
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assume that the b quark Yukawa interaction must be kinematically suppressed in 
order to avoid thermalization of the condensate on a time scale small compared with 
H _1 and so to allow the asymmetries to form. 
The condition < ip > > kT v corresponds to 

A 2 / TT \ 3/2 

T R <-^( K + m ,)M P1 )^(-) (5.3), 
where during matter domination, 

with a Q the scale factor at H m m s . Using equation (4.11) together with the observed 
B asymmetry, ^f- ks 10~ 10 , we find that 

R ~ ^2nsg5 p Kk r V m s / VmJ 



Thus, with Ay, = At, ~ 4x10 2 for the b quark Yukawa coupling and with k = 30, 
we find that the condensate will survive if 

T^f^fViLfVoeV (5.6), 



m, / \ m. 



where we have used m s 10 2 GeV. The form of this constraint depends on whether 
we impose the no-evaporation constraint before the asymmetry has fully formed 
or not. If we consider the constraint to apply at H « m s , when the asymme- 
try is minimal, then the T R upper bound is proportional to ^ • On the 
other hand, if we allow the ip asymmetry to grow to its maximum value before 

/ 111 fk -f- III \ 5 / 1 6 

thermalization, then the Tr upper bound is proportional to ( — ^ — -J . In both 
cases the upper bound on Tr is weakened by having m^ > m s . We see that, with 
iru < m s , T R can take values up to around lo7 gg v without preventing the formation 
of the asymmetry. With m^ > m s this constraint becomes weaker, allowing any 
reheating temperature up to the thermal gravitino limit to be compatible with the 
no-evaporation constraint, regardless of the value of 9, so long as m^ is sufficiently 
large. 
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We also have to check that the if) condensate does not decay before the conden- 

/ \l/4 

sates fully form at H pa [ m ^ m J ni s . If < if) > > m s , then the two-body 
tree-level decay will be kinematically suppressed. This occurs if 

> ( m s \ 1 / 2 



Vx ~ vm p1 ; (5 - 7) ' 

which will almost certainly be satisfied. The higher-order decay modes will then be 
suppressed by a factor of at least ( A <#> ) ' which gives r d C H at H k m s . Thus 
if) decay is ineffective at H pa m s and the no-evaporation constraint is the correct 
condition for the initial if) asymmetry to be able to form. 

A second, perhaps less important, constraint on the reheating temperature comes 
from the requirement that — ~ 1CT 10 can be consistent with non-renormalizible 
operators whose natural mass scale is the Planck scale, corresponding to |A| < 1. In 
fact, from (4.11), we see that, with ^ pa 1CT 10 , the reheating temperature is given 
by 

T R Pa J*L ( ^t±^lY /2 109Ge V (5.8). 

(The conventional ip A d=4 models give the same result but with m^ — > and 8 P ~ 1 
[§]). Thus we see that |A| ^ 1 is necessary in order for Tr to be consistent with 
the thermal gravitino bound. Although |A| ^ 1 is possible even if the mass scale of 
the non-renormalizible operators is small compared with Mpi, it is most natural for 
the case of Planck scale operators. Thus d=4 models are most naturally consistent 
with the thermal gravitino bound when the mass scale of the non-renormalizible 
operators corresponds to the Planck mass. For the case of a conventional single 
field A-D mechanism with d=6 operators, the reheating temperature is given by [|J 

T R pa ( U_ j lOGeV (5.9). 

Thus in this case the reheating temperature is expected to be very low compared 
with the thermal gravitino bound for the case of Planck scale non-renormalizible op- 
erators, although for non-renormalizible operators with a smaller mass scale, which 
would naturally have |A| ^> 1, larger reheating temperatures would be possible. 
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Therefore we see that d=4 models are favoured for the case of Planck scale oper- 
ators, naturally allowing for a much wider range of reheating temperatures than 
d=6 models, whilst d=6 models are favoured if the natural mass scale of the non- 
renormalizible operators is much smaller than the Planck scale. 
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6. Condensate Thermalization and Decay after Reheating and an 
Upper Limit on Dirac Neutrino Masses. 

After the inflaton decays, the Universe will be radiation dominated with H = 
where k H = ( ^45^ ) ^ ~ 17- We first show that tp will typically thermalize 
at a temperature large compared with the temperature of the electroweak phase 
transition T EW « 10 2 GeV. T s > H occurs if Ay, > (^f^ (m^) ^ Since 
T R < 10 10 GeV, this will be satisfied if Ay > 5xl(r 2 (^f /A ". Thus so long as the 
scattering process in not kinematically or Boltzmann suppressed, corresponding to 
A ' / '^ [ f > smaller than 1, the condensate will thermalize. With, for T < T R , 

< * >ss aw (loio^vJ l-5T"J UJ (fu) ' 

where we are using m s m 10 2 GeV throughout, A,/ '^ > < 1 requires that 
Ay, < in4 / m s \ 1/4 /!0 10 GeV\ 1/2 /mA 1 / 2 



This will typically be satisfied for some value of T larger than T E w- 

For the case of the gauge singlet (p condensate, we first note that since the <fi scalar 
is R p odd, it can decay only if it is not the lightest supersymmetric partner (LSP). 
The most rapid possible <f> decay will correspond to a tree-level two-body decay to a 
left-handed neutrino and a neutralino via the neutrino Yukawa coupling \ u . This will 
be kinematically allowed so long as one of the neutralinos has a mass less than the 
4> scalar mass. In particular, this will occur if one of the neutralinos is the LSP, as is 
strongly favoured by the possibility of neutralino cold dark matter. The (p decay rate 
will then depend on the proportion of light mass eigenstate neutralino(s) contained 
in the weak eigenstate fermion in H u . This will in turn depend on whether the 
light neutralino in question is mostly gaugino or Higgsino. If it is mostly Higgsino, 
then the decay will occur via the neutrino Yukawa coupling with coupling strength 
approximately equal to A„. On the other hand, if it is mostly gaugino, then we would 
expect the coupling of the scalar to the light neutralino to have an addtitional 
suppression factor of the order of where a corresponds to the yuH u H d term in 
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the MSSM superpotential [|IJ. Typically, for /x < ITeV, this factor will not be much 
smaller than about 0.1. 

Thus in the case of a mostly Higgsino light neutralino (or, more generally, for 
the case where the mass term \x is small compared with the scalar mass), tree- 
level two-body (f) decay will occur via the neutrino Yukawa coupling \ u if « 
(m^ + m;?) 1 / 2 > H, where a v « ^. This is satisfied if 

The decay is kinematically allowed so long as \ p < <fi > < (m| + m 2 ) 1 / 2 . With, for 



T < T 



K ■ 



,d>~ T R m 2 (m s 2 + m s m,)3/ 4 /M Pi y/ 2 /Tx3/ 2 



this condition requires that 



K < 1n4 /!0 10 GeV\ 1/2 /m^ + m s \ 13 / 4 /m s \ 3 / 2 



It is straightforward to show that, in the cases of most interest to us here, the decay 
of the cf) condensate typically occurs before it can thermalize by scattering. The 
condition for the condensate to thermalize by scattering via the Yukawa coupling 
\ v is that T s « k r aA 4 T > H = where a « |log (5^5:) ■ This implies that 

/ np x 1/4 / 1 x 1/4 

A,>5xl0- 4 ( — ) (-) (6.6). 



m s / \a 



This lower bound is typically larger than the lower bound on \ u coming from (p decay. 
It is possible that thermalization could occur by scattering from light sleptons and 
SU(2) gauginos in the plasma, which would replace X u by (A^g) 1 / 2 , where g « 0.6 
is the SU(2) gauge coupling. However, the lower bound on X u will still typically be 
larger than that coming from decay. In particular, this is true for the important 
case of (p decay below the temperature of the electroweak phase transition, which 
we discuss below. Another possibility for thermalizing the condensate is via inverse 
decays and related 2 — > 1 processes, which are expected to have a rate r inv « kA 2 T, 
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where k< 1. Although at high temperature this rate can be large compared with 
the decay rate, for temperatures at or below the electroweak phase transition 
temperature, which are of most interest to us here, the direct decay rate will be 
much larger than the rate of thermalization via inverse decays. In all this we have 
assumed that the decay occurs via the neutrino Yukawa coupling. It is also 
possible that the condensate could decay via the non-renormalizible superpotential 
coupling ^0u c d c d c once < ip > is introduced, which gives an effective 0-0-0 coupling. 
However, it is straightforward to check that this effective coupling is in general much 
smaller than the typical values of X u considered in neutrino mass models, and so 
may be neglected when discussing condensate decay. 

Thus we see that the condensate can evade decay until T < Tew if K is suffi- 
ciently small. To see what this implies for the baryon asymmetry and for neutrino 
masses, we first note that the gauge singlet scalar will typically correspond to a 
linear combination of the three right-hand sneutrino generations. Thus X u will typ- 
ically correspond to the largest neutrino Yukawa coupling up to a factor of around 
-^=. In the m^ — > limit, corresponding to the case where the neutrinos have Dirac 
masses, we can define an unbroken B-L asymmetry by defining to have B=l. How- 
ever, so long as the condensate decays after the electroweak phase transition has 
occured, the effect of anomalous electroweak B+L violation ||, which is in thermal 
equilibrium at temperatures larger than Tew 5 will be to alter only the B asymmetry 
coming from the thermalized -0 condensate and so prevent a cancellation of the B 
asymmetry coming from the if) and conensate, even though the net B-L asymme- 
try will be zero. Thus so long as the condensate decays at Td </, ^ Tew in the limit 
m</> — ► 0, there will still be a B asymmetry. The magnitude of the B asymmetry 
will be essentially the same as that previously calculated from the decay of the ip 
condensate alone, since the magnitude of the net B asymmetry will, up to a factor 
of the order of 1, equal that of the B-L asymmetry coming from ip thermalization at 
T Z Tew- The coupling X u is related to the heaviest neutrino mass in the m^ — > 
limit by m„ ps 10 2 A^GeV. From this and (6.3) we find that, for the case of a mostly 
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Higgsino light neutralino, the </> condensate decays at a temperature Td^ ~ I0 7 m u . 
Thus the condition T<j ^ Tew implies that all Dirac neutrino masses should satisfy 

m„ < lOkeV (6.7), 

assuming that </> corresponds to a roughly equal combination of the three sneutrino 
generations, as we would generally expect. This is true for the case of a light 
neutralino which is mostly Higgsino, or more generally for the case where the <fi 
decay via the Yukawa coupling X u is completely unsuppressed. From this we see that 
so long as the Dirac neutrino masses are all below about lOkeV, the condensate 
will generally decay below the temperature of the electroweak phase transition and 
so a baryon asymmetry will be generated even in the limit of B-L conservation. On 
the other hand, for the case of, for example, a mostly gaugino LSP this upper bound 
could be increased to around lOOkeV or more, depending on the particular gaugino 
LSP mass eigenstate and the \i parameter. These upper bounds should be compared 
with the present experimental upper bounds on the neutrino masses, m„ T < 24MeV, 
m. v < 160keV and < 5.1eV [I3|. From these we see that the requirement that 



a non-zero B asymmetry can be generated in the limit of unbroken B-L in the case 
where the LSP is a neutralino imposes a non-trivial upper bound on m- and m Ur . 
In particular, we see that it would be possible, in principle, to experimentally rule 
out this class of Affleck-Dine models, for example if neutrinos with Dirac masses 
significantly larger than around lOkeV were found to exist together with an LSP 
which was mostly Higgsino. We also note that an unbroken B-L asymmetry would 
rule out the possibility of the d=4 LH U direction, leaving the d=4 u c d c d c direction 
as the unique d=4 possibility in the case of B-L conserving models with non-zero 
neutrino masses. 

For the case with m^, ^ 0, the neutrinos will gain Majorana masses via the 
see-saw mechanism fli~4" |, with m u w ( A " 100GeV ) _ Thus in this case Td</> is given by 



Therefore typically the <fi condensate will decay at T > Tew i n this case. 
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Throughout the above discussion we have assumed that the Universe is radia- 
tion dominated. It is straightforward to show that this is indeed the case. The if) 
condensate would dominate the energy density only once T satisfies 

,2 



(in: nijii,;,) 1 2 'IV: 



(6.9), 



A Mpi 

which is typically satisfied only for temperatures less than around lCT 7 GeV. For 
the <fi condensate the energy density is even less than the ip condensate, by a factor 
(—^r — ) ■ Thus the Universe will be radiation dominated when the condensates 
thermalize or decay. 
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7. Conclusions 

We have considered the possibility of generating the observed baryon asymme- 
try via an Affleck-Dine mechanism based on the renormalizible F- and D-flat u c d c d c 
direction of the SUSY Standard Model. In order to avoid breaking R p whilst al- 
lowing a d=4 superpotential term to lift the flatness and drive baryogenesis, we 
considered extensions of the SUSY Standard Model which have additional gauge 
singlets <fi, such as commonly occur in models which seek to account for neutrino 
masses. In such models the u c d c d c direction becomes potentially as important as 
the more commonly considered LH U direction. We have shown that the A-D mecha- 
nism based on the d=4 operator (p^ 3 , where ifj is a gauge non-singlet A-D field, can 
indeed (for an appropriate choice of parameters) generate the baryon asymmetry 
whilst allowing the scalar fields to evolve to a phenomenologically acceptable min- 
imum. The resulting asymmetry is suppressed relative to the asymmetry coming 
from the more conventional ^ 4 -based A-D mechanism (such as the LH U direction) 
by a factor (^ m "^ s m J , ah couplings and CP phases being taken equal, where m^ is 
the SUSY (j) mass term and m s is the soft SUSY breaking mass scale. This suggests 
that m^ cannot be much larger than m s , if the observed baryon asymmetry is to be 
generated without requiring very small couplings in the non-renormalizible terms. 
The requirement that the initial condensate particle asymmetry can form before the 
condensate is thermalized imposes an upper bound on the reheating temperature of 
10 ^ cV in the limit where m^ < m s , where 9 is the CP violating phase responsible 
for the baryon asymmetry. The upper bound becomes weaker if m^ > m s . Thus 
with a sufficiently large m^ or small 9 the whole range of reheating temperatures up 
to the thermal gravitino constraint can be compatible with the initial formation of 
an asymmetry. The -0 condensate will typically thermalize before the electroweak 
phase transition occurs. Then in the limit of unbroken B-L (for which case the 
u c d c d c direction is the unique d=4 possibility), which corresponds to m^ — > 0, a B 
asymmetry can be generated only if the 4> condensate decays below the temperature 
of the electroweak phase transition, when anomalous electroweak B+L violation is 
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out of thermal equilibrium. This will generally be true if all neutrino masses are less 
than around lOkeV. In the case where the LSP is a neutralino, or more generally 
where there is a neutralino mass eigenstate lighter than the <fi scalar, the decay 
condition imposes a non-trivial upper limit on Dirac neutrino masses. For example, 
for the case of a mostly Higgsino LSP, the upper bound is around lOkeV, whilst 
for a mostly gaugino LSP this bound could increase to around lOOkeV or more, 
depending on the \x parameter of the MSSM and the particular gaugino LSP mass 
eigenstate. Thus the observation of a Dirac mass for the ji or r neutrino signifi- 
cantly larger than lOkeV together with a mostly Higgsino LSP, for example, would 
experimentally rule out this class of Affleck-Dine models. 

The author would like to thank Kari Enqvist for useful comments. This research 
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